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The enormous distances over which astrophysical jets propagate without los- 
ing their coherence certainly constitute one of the most striking features of 
these objects. Typically, jets from Young Stellar Objects (hereafter YSOs) 
do reach out to a few parsecs, while the radial extent of their region of ori- 
gin appears to be smaller than ~ 100 A.U, making jets extremely elongated 
structures. 

Blandford and Rees [6] already pointed out in their 1974 pioneering work 
that in laboratory experiments, jets do not propagate much farther than about 
ten times their radii, which makes the propagation lengths of astrophysical 
jets all the more impressive. The simplest way out of this conundrum would 
be to assume that jets are ballistic. Indeed, for YSO jets at least, the observed 
opening angle (~ 5°) is consistent with the idea that they freely expand when 
one compares their thermal and bulk velocities. However, this option leaves 
open the issue of the formation of such powerful jets in the first place. And, 
as critically, the ballistic hypothesis does not explain how these jets survive 
the development of the Kelvin-Helmholtz instability, which is now known to 
be quite disruptive in purely hydrodynamic jets [7] [8]. 

The shortcomings of the simple ballistic picture certainly motivated to 
some extent the elaboration of MHD jet models. Such models, however, are 
also prone to instabilities. The most important ones discussed in the literature 
can be grouped into three categories: 

• MHD Kelvin-Helmholtz instability. As for its HD counterpart, the driving 
agent is the velocity gradient at the jet/external medium interface. This 
instability has received a lot of attention in the literature, as the largest 
source of free energy in a jet is its bulk motion (see [5] and [H] for reviews). 

• Conversely, the presence of a magnetic field provides a source of instability 
even in the absence of bulk motion. Ideal MHD instabilities are commonly 
divided into current- and pressure-driven, according to the driving factor 
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(equilibrium parallel current in the first case, and gas pressure versus field- 
line curvature in the second). 
• Radiative instabilities, related to the coupling of the radiation field and of 
the plasma dynamical quantities. 

The structure of jets is not precisely known, which is one of the difficulties 
in analyzing their stability. Most stability analyzes assume that jets can be 
described as some sort of cylindrical column in motion, pervaded by a mag- 
netic field. Self-collimated jet models are not exactly cylindrical, but as the 
observed opening angles are small, the assumption of cylindrical shape is not 
expected to be a major limitation. 

More critically, such jet models have a helical field structure, with the az- 
imuthal component of the magnetic field dominating over the vertical one in 
the outer jet regions. This follows in most models because the magnetic ten- 
sion is the confining force ensuring self-collimation. Static MHD columns (i.e., 
not subject to the bulk motion characterizing MHD jets) pervaded by a helical 
magnetic field are referred to as "screw pinches" in the fusion literature. It 
is also known in this context that the dominance of the azimuthal field com- 
ponent leads to both types of MHD instabilities mentioned above, and may 
cause the disruption of the plasma column itself on a few dynamical time- 
scales. This has long been an argument against magnetically self-confined jet 
models. However, recent investigations indicate that a bulk motion can play 
an important stabilizing role (see section 5 for pressure-driven instabilities, 
and e.g. [1] and references therein for current-driven instabilities). Conversely, 
the presence of a magnetic field can help stabilizing the Kelvin-Helmholtz 
modes [53], [2]. These recent advances seem to indicate that a sophisticated 
equilibrium jet structure is required if one is to understand jet stability prop- 
erties, a state of development not yet reached by the subject, but that now 
appears to be within sight. 

To conclude these introductory remarks, I would like to point out that, 
in the nonlinear phase, an instability can have three broad types of outcome: 
i/ disruption of the fluid configuration (in the case at hand, of the jet as a 
jet); ii/ internal reorganization, the flow becoming laminar again in the end; 
iii/ turbulence (with or without internal reorganization of the structure). The 
most prominent objective of the study of jet stability is to understand how the 
first issue is avoided in real jets; this issue may well be seen as our inability 
to formulate the initial value problem correctly. A second but important issue 
is to understand how turbulence might be driven by jet destabilization. This 
issue is probably more important in AGN than in YSO jets, as turbulence 
is often invoked in the former context as a source of high energy particle 
acceleration. 

The object of these lecture notes is pressure-driven instabilities. As most 
investigations of this problem have been made in the fusion context for 
static columns, this essential aspect of the subject will first be reviewed, be- 
fore briefly presenting the more recent (and more scant) results on moving 
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columns. The next section presents some general ideas about the physical ori- 
gin of MHD instabilities; the concept of magnetic shear is introduced there 
as well, and its stabilizing role on pressure-driven instabilities, expressed by 
Suydam criterion, is discussed. Section [3] introduces the Lagrangian form of 
the perturbation equations used in static columns, as this formulation is the 
most powerful to derive general results, such as those derived from the Energy 
Principle presented in section [3] Section [4] presents the dispersion relation of 
pressure-driven instabilities in low magnetic shear that are expected to char- 
acterize jets outer regions. The few published results on moving columns (i.e., 
jets) are presented in section [5] The last section summarizes the present state 
of understanding of this aspect of jet stability, and outlines areas where im- 
provement is needed. 

Sections Q] and [5] are intended for a general audience, while sections [2j [3] 
and[5]are more theoretical in nature. The exposition is aimed at the graduate 
student level. 



1 Heuristic description of MHD instabilities 

This first section is intended to provide the reader with some qualitative and 
semi-quantitative ideas about the onset and characteristics of pressure-driven 
instabilities, leaving technical aspects of the stability analysis to the later 
sections. 

1.1 Qualitative conditions of ideal MHD instabilities in static 
equilibria 

The equilibrium configurations leading to an ideal MHD instability have been 
well investigated in the fusion literature. For current-driven instabilities, the 
first criterion was devised by Kruskal and Shafranov. Basically, it states that 
in cylindrical column of length L, instability follows if the magnetic field line 
rotates more than a certain number of times around the cylinder, from end to 
end. The exact number of rotations required for instability is dependent on 
the considered equilibrium configuration; it is usually of order unity. 

Concerning pressure-driven instabilities, a more clear-cut necessary con- 
dition of instability can be stated: instability follows once the pressure force 
pushes the plasma outwards from the inside of the field line curvature. This 
condition can be derived from the Energy Principle, as will be shown later on. 

These conditions of onset of instability are illustrated on Fig. [TJ In an 
actual plasma, the origin of an instability (current- or pressure-driven) is 
usually not easy to pinpoint except in special instances. For example, if the 
plasma is cold (no pressure force), the instability is necessarily current-driven. 
Also the growth rates of current-driven modes are known to decrease with 
spatial order - e.g., they decrease with increasing azimuthal wave-number m 
- while the most unstable pressure-driven modes have a growth rate which 
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Fig. 1. Qualitative description of the conditions of onset of MHD instabilities (see 
text for details). 



is nearly independent on the wavenumber. Consequently, large wavenumber 
unstable modes are therefore always pressure-driven in a static, ideal MHD 
column. Besides these two limiting cases, an MHD instability almost always 
results from an inseparable mix of pressure and current driving. If the column 
is moving, the distinction between Kelvin-Hclmholtz, current- and pressure- 
modes is even more blurred, except in some cases, where branches of instability 
can be identified by taking appropriate limits. 

In terms of outcome of the instability, it is essential to know whether 
unstable modes are internal or external, i.e., have vanishing or substantial 
displacement on the plasma surface (here, the jet surface) . It is well-known 
in the fusion context that unstable external modes are prone to disrupt the 
plasma, as may be the case, e.g., with the m = 1 ("kink") current-driven 
mode. 

In the next sections, mostly high wavenumber modes will be examined, 
where the pressure-driving is most obvious, in order to best identify the char- 
acteristic features of this type of instability. 



1.2 Magnetic shear, magnetic resonances, and Suydam's criterion 

The concept of magnetic shear plays an important role in the understanding 
of the stability of pressure-driven mode. The magnetic shear characterizes the 
change of orientation of field lines when moving perpendicularly to magnetic 
surfaces. In the case of cylindrical equilibria, this concept is illustrated on 
Fig. [21 Magnetic surfaces are cylindrical. Field lines within magnetic surfaces 
have an helix shape; the change of helix pitch rB z / Bg, characterizes the mag- 
netic shear. A quantity related to the pitch and largely used in the fusion 
community is the safety factor q: 

(i) 



R Bg 

where R a is the column radius. For reasons soon to be discussed, a high enough 
safety factor is required for stability, hence its name. The magnetic shear s is 
defined as 
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r dq 
q dr 



(2) 




Fig. 2. The change in the pitch of field lines between magnetic surfaces is the source 
of the magnetic shear (see text). 

Magnetic resonances constitute another important key to the question of 
stability. A cylindrically symmetric equilibrium is invariant in the vertical and 
azimuthal direction, so that perturbations from equilibrium can without loss 
of generality be expanded in Fourier terms in these directions and assumed to 
be proportional to expi(m8 + kz). Magnetic resonances are the (cylindrical) 
surfaces where the wave vector k = m/reg + ke z is perpendicular to the 
equilibrium magnetic field: 



where k\\ is the component of the wave vector parallel to the equilibrium 
magnetic field. The significance of these surfaces stems from the fact that in 
general, dispersion relations incorporate a stabilizing piece of the form Vj/c|, 
where Va is the Alfven speed. This term is responsible for the propagation 
of Alfven waves, and arises from the restoring force due to the magnetic ten- 
sion (see section @] for the precise meaning of these statements) . As such, it 
is always stabilizing. Obviously, this stabilization is minimal in the vicinity 
of a magnetic resonance for a given (m, k) mode, so that pressure-driven in- 
stabilities are preferentially triggered at magnetic resonances for any given 
mode. 

Note however that a large magnetic shear limits the role of magnetic reso- 
nances in the destabilization of the plasma. Indeed, defining the perpendicular 
wave number 




(3) 
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k x = ~(-B g -kBg), (4) 
B \r J 

and designating by r c the radial position of the magnetic resonance of the 
(to, k) mode, one finds that 



to first order in (r — r c )/r c in the vicinity of the magnetic resonance r c . This 
implies that Fa&|| will remain small either if s -C 1 (small shear) or if the 
magnetic field is mostly perpendicular (\B$\ <C \B Z \) or azimuthal (\B Z \ <C 
\Bg\) so that \B e B z \/Bl < 1. However, if the field is predominantly vertical, 
it is little curved, and pressure destabilization is expected to be weak or non- 
existent according to the description of the condition on instability depicted 
in Fig. ??; furthermore, s being a logarithmic derivative is usually of order 
unity. Therefore, in practice stabilization by magnetic tension will be reduced 
essentially when the field is mostly azimuthal. 

These features are embodied in Suydam criterion, which expresses a suffi- 
cient condition for instability: 

Bl , dP 



V + — > 0. (6) 
8(U r ar 

The converse of this statement is a necessary condition for stability. The origin 
of this criterion is briefly discussed in section[31 It turns out that this condition 
is both a necessary and sufficient condition of instability for large wavenumber 
modes [T3]. The condition of instability requires dP/dr < 0, which agrees with 
our heuristic description of the onset of instability given above. It will be also 
discussed in Section [3] that the growth rates 7 of pressure-driven instabilities 
are 7 ~ Cs/R (Cs is the sound speed and R a the jet radius). 

Coming back to Eq. ([6|), the first term is stabilizing, but the stabilization 
will be minimal in the condition just discussed, i.e., when the field is mostly 
azimuthal. Indeed, in this case, the equilibrium condition Eq. pip implies 
that dP/dr ~ Bgj\i r » B z /r ~ B z s 2 /r. This situation is expected to hold 
in magnetically self-confined jets outer regions. Indeed, most such jet models 
(e.g., [13] and [T5]) have \Bg\ 3> \B Z \ in the asymptotic jet regime to ensure 
confinement. This feature combined to the previous statement that MHD 
instabilities involving the boundary are most prone to disrupt static MHD 
columns makes the assessment of the role of pressure-driven instabilities in 
MHD jets particularly critical for the viability of such models. This viability 
hinges on the hopefully stabilizing role of the jet bulk motion (see section EJ. 
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2 Ideal MHD in static columns: 

The simplest framework in which the stability of jets can be investigated is 
ideal magnetohydrodynamics (MHD). Justifications and limitations of this 
approach are briefly discussed in Appendix 1X1 



2.1 Equations 

The MHD equations used in these notes are the continuity equation, the 
momentum equation without the viscous term, the induction equation without 
the resistive term, and a polytropic equation of state. Incompressibility is 
not assumed, as pressure-driven modes are not incompressible except at the 
marginal stability limit. These equations read 

^ + V P v = 0, (7) 

dv _ VP T B VB 

+ v .y. v = L + , ( 8 ) 

at p fi p 

3B 

— = Vx(vxB), (9) 

P = Kp\ (10) 

with standard notations, and where Pt = P + B 2 /2p, Q is the total (gas and 
magnetic) pressure, K a constant, and 7 the polytropic index. 



2.2 Equilibrium 

Using a cylindrical coordinate system (r,9,z), a static (v = 0) cylindrical 
column of axis z is described by a helical magnetic field Bg{r), B z (r), and a 
gas pressure P(r) depending only on the cylindrical radius r. The continuity 
and induction equations are then trivially satisfied, as well as the vertical and 
azimuthal component of the momentum equation, while the radial component 
reduces to 

dP T Bl , N 

-_Z-_L = 0. (11) 

dr p, r 

This cylindrical equilibrium is best characterized by introducing a number 
of quantities homogeneous to an inverse length, both in vectorial (K-b, JCp 
and Kc) or algebraic form (ICb, K-p and tCc)- They are defined by: 



(12) 
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VP 1 dP 
Kp = -p- = -p- Q ~d^ er = /Cpe '' (13 ) 

B 2 

K c = e,| • Ve|| = -^e r = K c e r . (14) 

where B Q and P D are the equilibrium distribution of magnetic field and gas 
pressure, and e | = B Q / B Q is the unit vector parallel to the magnetic field; Kc 
is the curvature vector of the magnetic field lines, and Kb characterizes the 
inverse of the spatial scale of variation of the magnetic field, while Kp char- 
acterizes the inverse scale of variation of the fluid pressure. The first identity 
in these relations is general, whereas the second one pertains to cylindrical 
equilibria only. 

It is also convenient to introduce the plasma f3 parameter: 

/3^^p (15) 

This parameter measures the relative importance of the gas and magnetic 
pressures. 

With these definitions, the jet force equilibrium relation reads 

~K P = (K C -K B ), (16) 

Both forms of the equilibrium relation, Eqs. (jTTJ) and (116|) . express the 
fact that the hoop stress due to the magnetic tension (Kc) balances the gas 
(Kp) and magnetic (Kb) pressure gradient to achieve equilibrium and confine 
the plasma in the column. Self-confinement is achieved in this way when the 
external pressure is negligible at the column boundary. 



2.3 Perturbations: 



We want to investigate the stability with respect to deviations from equi- 
librium. As the background equilibrium is static, the problem is most easily 
formulated and analyzed in Lagrangian form: indeed, in this case, all equa- 
tions but the momentum equation can be integrated with respect to time. 
To this effect, we introduce, for any fluid particle at position r Q in the ab- 
sence of perturbation, the displacement £(r Q , t) at time t from its unperturbed 
position, so that its actual position is given by 

r(r ,t) = r + £(r ,t). (17) 
The unperturbed position r Q is used to uniquely label all fluid elements. 
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Denoting by SX the (Lagrangian) variation during the displacement of any 
quantity X, the linearized (Eulerian) equation of continuity d t Sp = — V(p i>) 
integrates intcQ 

5p = -V(p„£). (18) 
Similarly, the linearized induction equation dtSB = V x [B a x v) leads to 

SB = V x (B x £). (19) 

From these results and the polytropic equation of state, the total pressure 
variation reads 

5P T = -€ • VP - 7 P V • | + . (20) 

For static equilibria, one can without loss of generality take a Fourier 
transform the linearized momentum equation with respect to time. For a 
given Fourier mode, one can write £(r, t) = £(r) exp iut, so that the linearized 
momentum equation becomes 

-p u> 2 £ = -V5P T + 8T = F{g), (21) 

where ST = (B .WB + B.WB )//j, represents the variation of the magnetic 
tension forco The last identity in Eq. (|2"Tj) defines the linear operator F, 
operating on £ through Eqs. (fl9| and (|20|) . 



3 The Energy Principle and its consequences: 

The linear operator F of Eq. (|2~lj) is self-adjoint, i.e., taking into account that 
F is real: 

J V -F($)d 3 r = J |-F(77)d 3 r. (22) 

A demonstration of this relation can be found, e.g., in Freidberg [TB] (cf p. 242 
and Appendix A of the book) . 

As a consequence of this property of F, an Energy Principle can be for- 
mulated. Defining 

5W(?,S) = ~[e-m)d 3 r, (23) 

1 In these expressions, the difference between the Eulerian and Lagrangian vari- 
ations has been ignored as they disappear to first order in the displacement £ 
in the final equations. For the same reason, no distinction is made between the 
derivative with respect to r or r . 

2 Within a factor p . 
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and 

= P\t\ 2 d 3 r, (24) 
and taking the scalar product of Eq. (f2Tj) with leads to 

The self-adjointness of F has two important consequences (Energy Prin- 
ciple): 

• u! 2 is also extremum with respect to a variation of £. 

• Stability follows if and only if SW > for all possible £. 

Ascertaining stability through the last statement is usually an impossible 
task. Instead, one usually makes use of the Energy Principle in a less ambi- 
tious manner: if one can find some displacement making SW < then one 
has a sufficient condition of instability (or, taking the converse statement, 
a necessary condition of stability). This is actually how Suydam criterion is 
demonstrated. First the expression of SW is simplified by taking advantage of 
the cylindrical geometry and focusing on marginal stability and incompress- 
ible displacements (as they make SW more easily negative; see below). Next, 
one chooses a particular form of displacement in the vicinity of the magnetic 
resonance of an (m, k) mode, and looks under which conditions this displace- 
ment makes SW negative; the condition turns out to be Suydam criterion for 
a well-chosen displacement. These computations are rather lengthy and the 
reader is referred to Freidberg's book [16] for details. 

A useful fornjfl of SW has been derived by Bernstein et al. [I] , which reads 
(see Freidberg [IB] , p. 259) 



SW = - I (f'r 



2P (K P ■ S±)(K C • €1) - J\\(V X x e,|) • Q ± ] , (26) 



where £j_ is the component of the displacement perpendicular to the unper- 
turbed field B, Q± = V x (£j^ x B ) is the perturbation in the magnetic field, 
K-c is the curvature vector of the magnetic field and JCp the inverse pressure 
length-scale vector defined earlier; Jn and en are the current and unit vector 
parallel to the magnetic field, respectively. The quantities TCp and K-c are 
defined in Eqs. (H]) and (flgjl. 

The first term describes the field line bending energy; it is the term respon- 
sible for the propagation of Alfven waves, through the restoring effect of the 
magnetic tension, which makes field lines acting somewhat like a rubber band. 



3 The boundary term is ignored, as it is not required in this discussion. 
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The second term is the energy in the field compression, while the third is the 
energy in the plasma compression. The fourth term arises from the perpendic- 
ular current (as VP = x B in a static equilibrium), and the last one arises 
from the parallel current Jm . Only these two terms can be negative, and give 
rise to an instability if they are large enough to make uj 2 < 0. Pressure-driven 
instabilities are driven by the first of these two terms, while current-driven 
instabilities are due to the second one. Pressure-driven instabilities are further 
subdivided into interchange and ballooning modes, depending on the shape of 
the perturbation, but the basic properties of these different modes are similar, 
and this distinction will not be discussed further in these notefl 

For our purposes here, we are mostly interested in what can be learned 
from the form of the fourth term. First note that this term is destabilizing in 
cylindrical geometry when K-cK-p > 0; this justifies the necessary condition of 
instability given in section ITTTI Furthermore, Eq. ([25)) and ([^l imply that the 
pressure-driving term produces an inverse growth rate 7, of order of magnitude 

| 7 | 2 ~ C 2 s KclCp ~ C 2 s /Rl (27) 

where C5 is the sound speed and R Q the column radius. This is quite fast, 
comparable to the Kelvin-Helmholtz growth rate in YSO jets. This order of 
magnitude will be used in the next section to set up an ordering leading to 
analytically tractable dispersion relations for pressure-driven unstable modes. 

4 Dispersion relation in the large azimuthal field limit: 

Most of the general results on pressure-driven instabilities were obtained in 
the fusion literature either from the use of the Energy Principle, or from the 
so-called Hain-Lust equation (a reduced perturbation equation for the radial 
displacement [19) [17j). These approaches are quite powerful, but not familiar 
to the astrophysics community, and involve a lot of prerequisite. 

It is more common in astrophysics to grasp the properties of an instability 
through the derivation of a dispersion relation. There are actually two papers 
doing this in the jet context for pressure-driven instabilities; however, the 
first one, by Begelman |3J, focuses on the relativistic regime which brings a 
lot of added complexity to the discussion, and the second one [21] is partially 
erroneous. 

Fortunately, in the limit of a near toroidal field of interest here, a dispersion 
relation can be derived ab initio by elementary means, and this approach is 
adopted here. To this effect, it is first useful to reexamine the behavior of 
the three MHD modes in an homogeneous medium, in the limit of quasi- 
perpendicular propagation. It is known that this limit allows the use of a 
kind of WKB type of approach in the study of interchange and ballooning 

4 In particular, Suydam criterion applies also to ballooning mode in cylindrical 
geometry; see Freidberg's book [T5], pp. 401-402 for details 
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pressure-driven modes (see, e.g., Dewar and Glasser [H]), a feature we shall 
take advantage of in these notes. 

4.1 MHD waves in quasi-perpendicular propagation in 
homogeneous media: 

We consider an homogeneous medium pervaded by a constant magnetic field 
B . The analysis of linear perturbations in such a setting leads to the well- 
known dispersion relation of the slow and fast magnetosonic modes and the 
Alfven mode. Our purpose here is to point out useful features of these modes 
when the wavevector is nearly perpendicular to the unperturbed magnetic 
field. 




Fig. 3. Definition of the reference frame (eii,e;,eji) 

To this effect, let us consider plane wave solutions to Eq. (|2Tj) . where 
£ oc exp(— ik ■ r), and assume that the direction of propagation is nearly 
perpendicular to the magnetic field, i.e. <C k± (defined in Eqs. [5] and 
[4}. The focus on quasi-perpendicular propagation comes from the remarks 
of section 11.21 where it was noted that instability is easier to achieve in the 
vicinity of magnetic resonances, i.e., where ku <C k±. 

Let us also introduce the orthogonal reference frame (eii, e;, &a) where 
e | = B /B is parallel to the unperturbed magnetic field, e; = kj_/k±, and 
eA = eii x ei (see Fig. [3]). With our definition of fcii and k± in Eqs. © and 
([4}, ba = e r - The subscripts I and A stand for longitudinal and alfvenic, 
respectively (en, ei, and eA are the directions of the displacement of purely 
slow, fast and alfvenic modes in the limit of nearly transverse propagation 
adopted here, as shown below). 

Denoting (£||,£z,£a) the components of the lagrangian displacement £ in 
this reference frame, the momentum equation Eq. (f2Tj) yields the following 
three component equations 
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(c 2 -C 2 fcf) £|| =Cffc|| k ± C, 



(28) 



(- 2 - eg fci 



vik 2 ) 6 = e|fc|, 



(29) 



(^ 2 - Vj A: 2 ) £ A = 0, 



(30) 



while the total pressure perturbation becomes 



SP T = ~i Po [(C 2 S + V 2 ) k±it + Cl fc||^||] . 



(31) 



Eq. (f3"0"|) gives the dispersion relation of Alfven waves, lo\ = Vjk 2 , 
which decouple from the two magnetosonic modes described by the remain- 
ing two equations. The solutions of the magnetosonic modes are easily de- 
rived and possess the following important properties. Characterizing quasi- 
perpendicular propagation with the small parameter e = |A|i/As_l| <C 1, these 
two equations imply w| ~ e§V|/((7j + V%)k 2 and ~ 0(e^|| ) for the slow 
magnetosonic wave, while oJp ~ (Cj + V^)k\ and £|| ~ O(e^i) for the fast 
magnetosonic one. 

Furthermore, the £j momentum component Ec(. (|29p combined with Eq. (|30|) 
and the ordering of the displacement component just pointed out implies that 
SPt — to leading order in e for the slow magnetosonic mode; note that 
the same property holds by construction for the Alfven mode. The cancella- 
tion of the total pressure for these two modes is essential from a technical 
point of view, and will lead to substantial simplification in the derivation of 
a dispersion relation performed in the next subsection. 

4.2 Dispersion relation and Kadomtsev criteria: 

Let us now come back to cylindrical inhomogeneous equilibria. Remember 
from section [3] that the pressure-driving term will contribute a destabilizing 
term cu 2 ~ C|/i? 2 to the dispersion relation. This term will be able over- 
come the stabilizing effect of the restoring forces of the Alfven and slow mag- 
netosonic modes only if VA|fc|||, Cs |fc||| < Cs/R . This constraint can be 
achieved in the vicinity of magnetic resonance as already previously noted. 

More precisely, a simplified dispersion relation can be found in the WKB 
limit with a displacement of the form £(r) = £ x exp —i(k r r + md + k z z), if 
the following ordering is satisfied: 

• |fc||r| <C or < 1 < \k r r\ <C \k±r\: the first inequality ensures that the 
stabilization by magnetic tension is ineffective (closeness to a resonance). 
The following inequalities ensure that a WKB limit can be taken. The 
implied ordering |A>ii| <C k± ensures that SP* will vanish to leading order 

5 For consistency with the previous sections, k± is the wavenumber in the longitu- 
dinal direction; it does not include the piece in the radial direction. 
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as in the homogeneous case discussed in the previous section. The last 
inequality allows us to neglect the contribution of the radial gradient of 
total pressure (which does not vanish), and greatly simplifies the analysis. 

• \B Z / Bg\ 2 s 2 \k± \ <C \ku |: this limit, which applies when \Bg | ^> \B Z \, ensures 
that the magnetic shear is not stabilizing. 

• | a; 2 1 <C |^f| 2 : this excludes the fast mode from the problem in the near 
perpendicular propagation regime considered here. As the fast mode is not 
expected to be destabilized in this regime (as \u 2 F \ 3> Vj/r 2 ), this does 
not limit the generality of the results while simplifying the analysis. 

It turns out that the resulting dispersion relation captures most of the 
physics of pressure-driven instabilities; this follows because the most unstable 
modes have growth rates nearly independent of the azimuthal wavenumber m 
|13j . and because current-driven instabilities are efficient only at low m and 
disappear from a WKB analysis. 

As previously, the projection Eq. (f2"Tj) on the longitudinal direction e; 
shows that the total pressure perturbation vanishes and that ~ |fc||/^J-l£|| ^ 
£|i ), while the components in the other two directions (en, e r ) are now coupled 
and read (some details of the derivation of these equations can be found in 
Appendix [BJ 

(u 2 - V 2 M kf) = - J -^1-V A 2 /Ccfc|t^, (32) 

where (3* = C 2 s /Vl, and V§ M = C 2 s Vl/{C 2 s + V%) is the slow mode speed 
in the near perpendicular propagation limit. The coupling of the modes blurs 
their character except in limiting cases. 
The quantity k 2 is defined as 

48* 

k ° = TTJ* c ~ 2/3 * /Cc/c "' (34) 

Note that if ICc = (i.e., when reverting to an homogeneous medium), 
Eqs. (f3"2"| and (f3"3"|) yield back the slow and Alfven mode, respectively. The field 
curvature couples the two modes. The quantity k 2 can be either positive or 
negative; the first term in Eq. (|34[) comes from the plasma compression, and 
the second one is the contribution of the pressure destabilizing term identified 
in section [3] 

As usual, these equations possess a non-trivial solution if their determinant 
is non zero, which yields the following dispersion for lo 2 : 



= *T J rT* y A Kck \\k> (33) 



(v; 



v£ M )kj + v 2 k 2 



viv 2 M 



k 2 (k 2 -28*ICclC p ) = 0. (35) 



Pressure-driven instabilities in astro-physical jets 



15 



First note that if both B z — (the so-called Z-pinch configurations) and 
m = 0, this equation is degenerate: one of the roots is uj 2 — and the other 
root is uj 2 = V\k 2 . Instability then requires that k 2 < 0, as k^ = in this 
case. This constrain is identical to the criterion^) derived by Kadomtsev from 
the Energy Principle for the m = mode in Z pinches (see Freidberg [16] p. 
286). 

When m ^ 0, Eq. ([35]) can be solved exactly but it is more instructive 
to analyze its properties. As the coefficient of uj 2 is equal to the sum of the 
two roots, and the last term is equal to their product, one finds that if k?, > 
2P*JCclC p , the two roots are stable, and if fcjj < 2/3*/Cc^p, one °f the roots 
is unstable. If B z = (Z pinch), this condition is identical to the criteriorfl 
derived by Kadomtsev for m ^ modes (see Freidberg jTC] pp. 284-285). 

Note that all these conditions of instability require ICclCp > 0, in agree- 
ment with the discussion of sections [3] and [TJ this condition is unavoidable in 
magnetically self-confined jets. The analysis presented here also shows that 
once this condition is satisfied, instability necessarily follows in static columns 
where \Bg\ ^> \B Z \ on some of the radial range, as the magnetic tension stabi- 
lizing effect V\k 2 is arbitrarily small in the vicinity of a magnetic resonance. 

Finally, the reader may ask how the local analysis presented here informs 
us on the global stability properties of the column. The answer lies in in 
oscillation theorem of Goedbloed and Sakanaka [18] . The theorem states that 
for any (m, k) unstable mode, the growth rate decreases when increasing the 
number of radial nodes. This implies that if an unstable mode with a large 
number of radial nodes is found (such as the modes considered here), an 
unstable nodeless mode will also exist, and this mode will have the largest 
growth rate. Such a mode will have a very disruptive effect on the plasma if 
its displacement is not vanishing on the boundary, as will be the case if the 
azimuthal field is dominant on the boundary. 

5 Moving columns: 

The previous section has shown that cylindrical columns with a predominant 
azimuthal magnetic field at least in some radial range are subject to pressure- 
driven instabilities. This situation holds in the outer region of self-confined 
magnetic jets, leading to a potentially disruptive configuration. However, in 
these regions a gradient of axial velocity due to the interaction of the moving 
jet with the outside medium is also expected to be present, and it is legitimate 
to investigate the effect of such a velocity gradient on the stability properties 
of pressure-driven modes. 



6 Kadomtsev's criterion for the m = mode in a Z pinch is a necessary and suffi- 
cient condition of instability, whereas the analysis presented here shows only the 
sufficiency of this condition. 

7 Same comment as in the previous footnote. 
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This problem has not yet been addressed in the astrophysics literature, but 
some relevant results are available in the fusion literature. In all the investi- 
gations cited below, the adopted velocity profile contains no inflexion point, 
in order to avoid the triggering of the Kelvin-Helmhotz instability. 

It is first useful to consider what becomes of Suydam criterion in presence 
of background motionqj. This investigation has been performed by Bondeson 
et al. [9]. Focusing on axial flows (U = U z (r)e r ), they conclude that the 
behavior of localized modes depends on the magnitude of 

M^ I/2 ^r, (36) 
<1 B z /q 

where the prime denotes radial derivative, and q is the safety factor (see 
section [2]). This quantity is a form of Alfvenic Mach number based on the 
velocity and magnetic shear, hence its name. When M 2 < f3, the flow shear 
destabilizes resonant modes. Above this limit, these modes arc stable, but in 
this case, unstable modes exist at the edge of the slow continuum, and may 
be global. The authors found however that in this case the growth rates are 
small (comparable to the resistive instabilities growth rates). Note also that, 
as q'/q ~ 1/r, M ~ {B$/ B z ){r / d)(U z /V A ) > 1 in MHD jets (d is the width 
of the velocity layer). 

These results seem to suggest that the region where the velocity shear layer 
takes place at the jet boundary is substantially stabilized in MHD jets. This 
seems to be confirmed by global linear stability analyzes, both for interchange 
and ballooning modes, except possibly for the m = ("saussage") mode [TU] 
[2"T] [25] . In all cases, increasing the flow Mach number efficiently reduces the 
amplitude of the displacement of the unstable modes at the plasma boundary, 
an important feature to avoid the disruption of the plasma. 

An efficient stabilization mechanism has also been identified in the non- 
linear regime by Hassam [20] . This author exploits an analogy between the 
m = pressure-driven interchange mode and the Rayleigh- Taylor instability 
in an appropriately chosen magnetized plasma configuration. From this analy- 
sis, he concludes that the m = pressure-driven mode is nonlinearly stabilized 
by a smooth velocity shear (dll z /dr ~ U/R ) if M s = U z /C s > [ln^/rg)] 1 / 2 , 
where r g is the instability growth time-scale (r g ~ c s (/C p /Cc') 1 ^ 2 ) and Td the 
diffusion time-scale (t^ ~ vK, p lCc where v is the viscosity, assumed compara- 
ble to the resistivity). The nonlinear evolution of an unstable, slightly viscous 
and resistive Z-pinch (i.e., a configuration where the field is purely azimuthal), 
was simulated by Desouza-Machado et al. [11]. They found that the plasma 
relaminarizes over almost all its volume for applied an velocity shear in good 
agreement with this analytic estimate. The core of the plasma still has some 
residual unstable "wobble" , which can apparently be stabilized by the mag- 
netic shear if some longitudinal field B z is added to the configuration. Note 



This requires a generalization of Eq. (|21|l ; also, the Energy Principle no longer 
applies as the resulting operator is not self- adjoint. 
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that the large values of Td/r g relevant to astrophysical jets lead to only weak 
constraints^ on the Mach number M s , so that this nonlinear stabilization 
mechanism is expected to be efficient in astrophysical jets. 

6 Summary and open issues: 

Pressure-driven instabilities occur in static columns when the pressure force 
pushes the plasma out from the inside of the magnetic field lines curvature, 
as shown from direct inspection of the "potential energy" of the linearized 
displacement equation (section [3]), and from a the dispersion relation of local 
modes (section 2]) . When unstable modes exist, the growth rates are of the 
order Cs/R where C$ is the sound speed and R a the jet radius. These are 
very large, comparable to the Kelvin- Helmholtz growth rate (the most studied 
instability in jets), especially that the ratio of the magnetic energy to the gas 
internal energy is expected to be of order unity (within an order of magni- 
tude or so). Such instabilities are known to be disruptive in the fusion context 
when the eigenmodes exhibit a substantial displacement of the plasma outer 
boundary; such a situation is relevant to magnetically self-confined jets, as 
the magnetic field in their outer region is predominantly azimuthal, a con- 
figuration most favorable to the onset of the instability (sections [T] and 2]) . 
However, the presence of a velocity gradient in the outer boundary due to the 
jet bulk motion is expected to have a substantial stabilizing influence, both 
in the linear and nonlinear regimes (section [5]) . 

In its present state, this picture possesses a number of loose ends: 

• The stabilizing role of an axial velocity gradient needs to be better under- 
stood. Not all modes may be stabilized in the linear regime, depending on 
the details of the equilibrium jet configuration, and the nonlinear mecha- 
nism identified in the literature is highly idealized and may not be generic. 
The one and only simulation of nonlinear stabilization published to date 
exhibits a very violent relaxation transient, which may still lead to jet dis- 
ruption. On the other hand, this transient is also an indication that the 
initial configuration of the simulation is way out of equilibrium, a situation 
which may not occur in real jets. 

• The role of jet rotation has not yet been correctly investigated. Preliminary 
results seem to indicate that it is stabilizing [21] ; however, jet rotation may 
not be an important dynamical factor in the asymptotic jet propagation 
regime. 

• Most investigations of pressure-driven instabilities rely on a very simple 
prescription of the equation of state, which raises an issue of principle. In- 
deed, the very large growth rates usually found for the instability indicate 
that it develops on time-scales much shorter than the collisional time-scale, 

9 For example Td/r g = 10 30 translates into M a > 8 only; in YSO jets, this ratio is 
most probably significantly smaller, and the constraint even weaker. 
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and the use of ideal MHD as well as a polytropic equation of state may be 
questioned in such a context, an issue briefly addressed in Appendix A 
more complex description of the plasma is required to validate the results 
obtained so far. 



Appendices 

A On the use of ideal MHD: 

In astrophysics in general, and jet stability analyzes in particular, an MHD 
framework is almost always adopted instead of the more precise kinetic one, 
due to its relative simplicity. The MHD approximation can be applied when 
the fluid is locally neutral, when all species can be described by a single fluid 
equation (i.e. when the relative drift velocity of species with respect to one 
another is small), and when Ohm's law is valid. The validity of these approx- 
imations has been discussed elsewhere [35] and will not be reproduced 
here; the interested reader is referred to these books for details. 

Furthermore, MHD stability (and jet stability in particular) is often inves- 
tigated within the framework of ideal MHD. Indeed, the dynamical time-scales 
of interest (including those of the considered instability) is almost always sub- 
stantially larger than the particle collision time-scale. Moreover, an isotropic 
pressure is often assumed, e.g. through a barotropic equation of state, and this 
raises another issue of principle within the framework of ideal MHD, as, in- 
deed, an isotropic pressure would be expected only if collisions at the particle 
level are not neglected. 

The isotropic pressure assumption can be justified to some extent by the 
fact that plasma microturbulence does limit pressure anisotropy to a factor of 
order unity. For example, within the framework of collisionless MHD, pressure 
anisotropy is self-limiting (for a recent synthetic discussion of this problem 
within the framework of collisionless magnetorotational instability, see Sharma 
et al. [25] and references therein). Nevertheless, this provides little support (if 
any) to the adoption of a closure in the form of a barotropic or adiabatic 
equation of state in a collisionless setting. 

Collisionless MHD approximations apply when the length-scales and fre- 
quencies under consideration are larger than the ion Larmor radius, and 
smaller than the ion cyclotron frequency, respectively. These conditions should 
be satisfied in jets, but I am not aware of any investigation of pressure-driven 
instabilities in this framework. Freidberg |16] argues that a simple rule of 
thumb to estimate the effects of the assumed closure is to replace the adia- 
batic index by 0, and to assume incompressibility of the motions within the 
framework of standard ideal MHD. 
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B Derivation of the dispersion relation: 

Some intermediate steps in the derivation of the dispersion relation of section 
H]are given here. The notations are the same as in this section. 

Direct computation of the pressure and magnetic field perturbation gives 

SP = -p C 2 s [V • i + IC P Cr} , (37) 

SB = B [V • (£ P + d) + (, r (K B + Kc)\ e,| 

-ik\\B {i r + ii). (38) 

This allows us to obtain the perturbation in total pressure and in magnetic 
tension: 

5P T = - Po {V\ + CD V ■ (£ r + €,) + ipoC 2 s k^\\ - 

Po^ ■ {VlK B + VlK c + C 2 S K P ), (39) 



ST = -VlK c [2V • (C P + €0 + 2£r ■ {Kb + Kg)] + 

ik\\V% [V • (€ r + €,) + 2^ r • /C c ] e ( | - fc|Vj(e r + (40) 

Furthermore, the equilibrium relation Eq. (jTSJ) allows us to eliminate K-b 
in terms of JCc and /C p = 7/Cp. 

The longitudinal component of the linearized momentum equation reduces 
to 8Pt = 0, once contributions of order fcy£z or £;/r are neglected in front of 
k±£_i. This constraint shows that £j ~ 0(fci|/fcxCr; fc||/fcxC||)- It a l so allows 
us to eliminate V • (£ r + from the remaining two component equations, 
which then reduce to Eqs. (|32|) and (|33|) . In the process, the contribution 
of dSPr/dr is shown to be negligible from the assumed ordering relations 
\k r \ <C and \B Z /Bg\ 2 s 2 \k±\ <C |fc|||, i.e., the magnetic shear stabilizing 
term can be neglected in this limit. 
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